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We study the interplay between a novel vortex-loop unbinding in finite magnetic field at T = Ty 
and flux-line lattice (FLL) melting at T = Tm in type-II superconductors. The FLL melts due to 
nucleation of vortex loops || c-axis, connected to flux lines. For moderate anisotropy, phase coherence 
1 1 c is lost at Ty > Tm due to an afe-plane vortex- loop unbinding with loops located close to thermal 
FLL fluctuations. For large anisotropy, phase coherence ||c is lost at TV < Tm due to nucleation of 
ab-plane vortex-loops uncorrelated to flux lines. 
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It is commonly assumed that the statistical mechanics 
of the flux-line lattice (FLL) in type-II superconductors is 
mainly governed by fluctuations of the flux lines around 
their positions in the ground state Abrikosov FLL, at 
least at magnetic fields well below the upper critical field 
such that vortex cores are well separated. Considerable 
progress in our understanding of the FLL has been made 
over the last few years, under this assumptiontl. 

However, we show that for large mass-anisotropies the 
low-energy thermal excitations in a superconductor in a 
low magnetic field are quite different. There exists a set 
of important topological excitations, responsible for de- 
stroying longitudinal superconducting phase coherence in 
a type-II superconductor, which are qualitatively differ- 
ent from fluctuations of flux lines. We also discuss in 
detail the fluctuations responsible for melting the FLL 
with no pinning. These two distinct types of fluctuations 
may influence each other in an essential way. This has ob- 
servable experimental consequences, as will be discussed. 
The induction range considered is < B/H C 2 <C 1, which 
is known to be the relevant limit for considering the phe- 
nomenon of FLL melting in extceme type-II superconduc- 
tors with uniaxial anisotropyfcm. Amplitude fluctuations 
of the superconducting order parameter are unimportant, 
and the anisotropic London model is appropriateu. 

The model we consider is a discretized, uniaxial 
anisotropic superconductor in the London limittj, on a 
lattice with L\ x L z in units of the numerical lattice 
constant d. The anisotropy is along the crystal c-axis, 
described by the parameter V = A c /A a = y/M z /M, 
where A a and A c are the magnetic penetration lengths 
along the crystal a6-plan and the crystal c-axis, respec- 
tively, and M z , M are normal-state quasiparticle effective 
masses along the z-direction and a6-planes, respectively. 
We take our coordinate system (x,y, z)-axes parallel to 
the crystal a, b and c-axes, respectively. Periodic bound- 
ary conditions Wi-iall directions are used. The Hamilto- 
nian is given byoQ 



h = 2tt 2 j J2 E ^-^KfcKfo), (I) 

where n^r,) is the integer vorticity through plaquette fj, 
at site Yi . Note that this model allows for arbitrary flux- 
line shapes, and also incorporates topological vortex loop 
excitations. Jo = $q(2/ '16ir 3 ; A 2 and Jo/r 2 are the energy 
scales for excitations involving n z (r) and (n x (r),n y (r)), 
respectively, and <I>o is the flux quantum. G(r) is the 
lattice London interaction, with Fourier transform 

r- 2 

Gx ' y( - k) = Ql + Ql + T-HQl+(PXa 2 ) 
r ^ _ (Q 2 + rf 2 A- 2 )/(Q 2 + d 2 A- 2 ) 
" W Ql + Ql + T- 2 {Ql+<P\- a 2 y U 

where — 2sin(fc M d/2), Q 2 = Q 2 . In these simula- 
tions, we have chosen A a = 0.75a, where a is the average 
distance between flux lines. 

To probe superconducting phase cohetence, we con- 
sider the helicity modulus T Al (g l/ ) (/i ^= vp. For a given 
perturbation in the external vector potential SA e ^ :t (q l/ )jl, 
"^((/j/) is the linear response coefficient giving the in- 
duced supercurrent j 

3M = -r it (g v ,T)8A^*(q v ). (3) 

We find that the helicity modulus in a uniaxial 
anisotropic lattice superconductor is given by 

T ll {q v ,T) _ _ 47T 2 J A 2 < n CT (q i ,)n CT (-g^) > 
T M (g„0) ~ VT 1 + (I + ^, 2 (r 2 -I))A 2 Q 2 ' 

(4) 

where < .. > denotes a thermal average, (n,v,a) are 
cyclic permutation of (x,y,z), and V is the volume of 
the lattice. 

We investigate the melting of the FLL W considering 
the structure factor for n z vortex elementsH 
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5(k) = <IS^ Z (r,)exp [ik-rj | 2 > (g) 

where iV z = X)i n z( r O- I n the ground state the FLL has 
the well know hexagonal form when B||c, and S(K, k z = 
0) has (5-function Bragg peaks at the reciprocal lattice 
points K. When the FLL melts, the Bragg peaks are 
washed out. The lowest T where S(K) vanishes thus 
defines the melting temperature Tm- 

We employ the following Monte Carlo (MC) procedure: 
We start with fixed average density f=B/$o of straight 
vortex lines parallel to the £-axis. We update the system, 
heating from the ground state by sweeping over the lat- 
tice in a systematic way. For each site we try to add one 
elementary closed vortex loop, choosing one of six pos- 
sible orientations of the vortex loop at random. These 
moves are accepted or rejected according to the standard 
Metropolis algorithm. This MC procedure provides a 
complete sampling of the phase space of the vortex vari- 
ables n(r), subject to two constraints: 1) A • B(r) = 0, 
where B(r) is the local induction and A is the lattice gra- 
dient operator. 2) The average induction is determined 
by B = (&o/V) J2i n ( r i) = f$oZ is constant. Note that 
our MC procedure allows for excitations of free vortex 
loops, i.e vortex loops not connected to any existing vor- 
tex lines. Each data point is obtained after discarding 
the first 5000 steps to allow for equilibration, whereas 
the subsequent 10000-40000 steps are used to obtain av- 
erages, where each step refers to a sweep through the 
entire L\ x L z lattice. 

In these simulations, we use / = 1/48, Lj_ = 24 and 
L z = 12. The reciprocal lattice vector of the vortex lat- 
tice K is taken to be the lowest in ^-direction consistent 
with the filling fraction / = 1/48 and a square numerical 
mesh, i.e. K = (4, 0)27r/L^. The q-value chosen in Eq. 
4 is the lowest finite one on our numerical lattice, i.e. 
q x = (1,0) 2tt/L ± and q y = (0,1) 2n/L ± . 

t x ,y will always be zero in finite field when the FLL 
is unpinned, since the FLL will move in response to an 
arbitrarily weak applied current. In our lattice model, 
the numerical lattice introduces artificial pinning which 
we must eliminate in order to be able to study the true 
melting of the FLL. This is achieved for low enough val- 
ues of / < f c , where we have found f c w 1/32 for the 
isotropic systems and decreasing very slightly with in- 
creasing mass-anisotropy. It suffices that T XiV vanishes at 
temperatures below those of which 5(K) vanishes. This 
indicates that the FLL 'floats' freely on the numerical 
lattice, having thermally depinned from it, in the form of 
an intact FLL as demonstrated by a finite S(K). This al- 
lows a study o£,the actual FLL-melting without artificial 
pinning effectsB. As shown in both Figs. |l| and ||, this is 
clearly so in our simulations. We have also checked that 
the T-dependence of T z , practically is not influenced by 
/ if / < 1/16. _ ' 

The quantitity Y z (q) essentially measures phase- 
coherence along the z-axis. Since vortex loops need not 
be connected to flux lines, T z (q) may vanish without 



a corresponding vanishing of S'(K). Thermal fluctua- 
tions may destroy superconducting phase coherence while 
leaving the FLL intact. This could possibly be inter- 
preted as the finite-field counterpart of the well known 
'inverted XY-transition' found in a_3-D lattice supercon- 
ductor model in zero magnetic fielcu. 
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FIG. 1. T,(g B ,T)/T,(«»,0) (A), T.(ft,,T)/T«( & ,0) (O), 
and 5*(K) (□) as a function of temperature for r = 1. Su- 
perconducting phase-coherence along the direction of the flux 
lines is left intact in the flux-line liquid S(K) = 0, indicating a 
melting transition into a disentangled vortex liquid. Shown in 
the inset is a snapshot of the flux- line system at T* < Tm- All 
excitations are thermally nucleated flux-line defects. There 
are no free vortex loops in the system, as explained in text. 
Note how T x vanishes before S(K); the system is in a 'floating 
solid phase'. 

In Fig. |l|, S(K), T x and T z are shown as functions of 
the temperature for T — 1. The hclicity modulus T 2 is 
seen to vanish at considerably higher temperatures than 
S'(K). Moreover, T x vanishes at a temperature distinctly 
below Tm , such that by the time we reach the FLL melt- 
ing temperature T = Tm, we have established a 'floating 
solid phase' from which the FLL can melt with no arti- 
ficial pinning effects to the numerical lattice. The fact 
that T z stays finite through the melting transition where 
S(K) vanishes, indicates that the flux lines stay intact in 
the vortex liquid, without large defects as evidenced by 
the still finite T z . S(K) vanishes due to loss of shear stiff- 
ness between flux lines with large effective tilt-moduli. 
As the temperature is increased further, T z eventually 
vanishes also, at T = Ty. For the isotropic case, this is 
primarily due to fluctuations in the flux lines, or equiva- 
lently directed vortex loops ||z-axis attached to flux lines. 
This is easily understood, since vortex loops ||l-axis and 
vortex loops in the afe-plane have the same self-energy 
when r = 1, but vortex loops ||z-axis attached to flux 
lines gain energy by annihilating one segment of the loop 
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with one segment of a flux line. Moreover, when vortex 
loops parallell to the a6-planes eventually are nucleated 
for r < 3, they tend to be nucleated in close vicinity 
to flux-line defects; they may interact attractively with 
the non-vertical (here horizontal) components of these 
defects. 
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FIG. 2. T z (q x ,T)/T z (q x ,0) (A), T x (q y ,T)/T x (q y ,0) (O), 
and S(K) (□) as a function of temperature for V = 16. Su- 
perconducting phase-coherence along the direction of the flux 
lines is lost before the FLL melts. Shown in the inset is a 
snapshot of the flux- line system at T* < Tm- Note the pro- 
liferation of afe-plane vortex loops uncorrelated to flux lines; 
the FLL is left intact. 

When T > 5, vortex loops in the a6-planes are easily 
nucleated anywhere in the a6-plane, due to the large re- 
duction of tfi£_sclf-cnergy of flux lines per unit length, by 
a factor 1 /rnlij. The result is a dramatic enhancement of 
vortex-loop unbinding and a resulting reduction of Ty at 
which T z (q) vanishes. This is shown in Fig. |^, where 
it is seen that, for T = 16, T z vanishes at temperatures 
well below the temperature Tm where S'(K) vanishes. T z 
vanishes, not due to fluctuations of the flux lines along 
the z-axis, but due to nucleation of vortex loops in the ab- 
planes located between flux lines. We conjecture that this 
transition very likely may be a finite-field counterpart of 
the 'inverted XY transition' first discussed by Dasgupta 
and Halperin of a lattice superconductor in zero-fielau. 
Flux-line defects are nucleated only at much higher tem- 
peratures, and the FLL is left completely intact at the 
vortex-loop unbinding transition, as shown by the value 
of S(K) when T z vanishes. This is also clearly illustrated 
in the inset of Fig. 0. 

Vortex-loops || z-axis also become easier to nucleate 
when r increases, but only by half the amount of the 
afe-plane loops, since an elementary loop || z-axis has two, 
and not four, flux-line segments parallell to the a6-plane. 
Moreover, it is always energetically favourable to attach 



them to already existing flux lines, by the same argu- 
ment as in the isotropic case. Hence, flux-line defects 
are also easier to nucleate in the anisotropic case, as one 
would have guessed from the reduction of-, the nonlocal 
tilt-modulus of the FLL when T increases^. Practically 
no vortex loops || z-axis are nucleated away from existing 
flux lines for T < Tm- 

For sufficiently large T, T z vanishes at temperatures 
well below those at which S(K) vanishes, i.e. Ty <C Tm- 
Nonetheless, the melting of the FLL is not necessar- 
ily 2D. Tm(T) continues to drop for increasing values 
of r which are much larger than the value of T where 
Ty = Tm- Thus, even if Ty <C Tm, Tm still depends on 
the coupling between planes. As long as this is the case, 
the FLL melting transition is therefore 3D in character. 
This is shown in Fig. 3, where both Ty and Tm are plot- 
ted versus T. Ty = Tm for T w 5. However, Tm oaly 
saturates to a finite value (which implies 2D meltingt3) 
at a much larger value, T ss 16. We have found an inter- 
mediate range of T-values where the entropy jump AS" 
associated with the FLL melting increases with Ftj. 
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FIG. 3. Melting temperature Tm (□) and vortex-loop un- 
binding temperature TV (A) as a function of T. Note satura- 
tion at r > 16, showing that the melting eventually becomes 
2D, but only for considerably larger values of F than those 
values where Tm ~ Ty- 

Previously, an entropy jump of AS*^ 0.3fcs per vor- 
tex per layer was found theoretically Li well in agree- 
ment with a number of experimentsEStJ. We now pro- 
pose a picture for the variation AS with T and induction 
B, when the perfect Abrikosov FLL has, been gauged 
away by a singular gauge-transformatioruJ. Hence, only 
thermally nucleated vortex-loops and thermally nucleated 
flux-line defects are relevant degrees of freedom, precisely 
as in our Monte Carlo simulations. 

The character of the vortex-loop unbinding and FLL 
melting transitions is governed by the effective interac- 
tions between vortex-loops in the afe-planes and the inter- 
action between flux-line defects, respectively. The effec- 
tive interaction between vortex loops is screened by ther- 
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mally nucleated flux-line defects, but is left unscreened 
by the flux lines in the rigid Abrikosov vortex latticSt3. 
Conversely, flux-line defects are screened by vortex loops. 
Note that this is different from the electromagnetic 
screening of the interaction between field-induced vor- 
tices, even in the perfect Abrikosov FLL. Such screening 
is present at all temperatures. 

Vortex-loop unbinding transition^ An unscreened in- 
teraction between vortex loopsiS makes the vortex- 
loop unbinding transition continuous, while screening 
will tend to drive it first order. Increasing T will en- 
hance thermal nucleation of flux-line defects, as argued 
above, enhancing screening of the vortex-loop interac- 
tion, driving the vortex-loop unbinding more first order. 
If one increases B, an increase in the density of flux- line 
defects will also enhance screening of the effective inter- 
action between vortex loops, enhancing the first order 
character of the transition. However, it will also sup- 
press fluctuations in individual flux lines, by increasing 
the flux-line tilt modulus C44. At low B C44 ~ B 2 , while 
the areal density of flux-lines n ~ B. Therefore, the 
dominant effect is the increase of the areal density of 
flux-line defects. Screening of the vortex-loop interac- 
tion is increased, driving the vortex-loop unbinding more 
discontinuous. There will be a crossover field B* above 
which the increase in C44 dominates, which will then re- 
duce the number of flux-line defects. Above this field, the 
vortex-loop unbinding transition will become less first or- 
der. 

FLL melting transition: An unscreened interaction be- 
tween flux-line defects will make the FLL melting tran- 
sition continuous. Screening of this interaction will tend 
to drive it first order. Increasing T will increase the 
nucleation of vortex loops in the a6-plane, enhancing 
the screening of the interaction between flux-line defects, 
driving the melting of the FLL more first order. On the 
other hand, if the magnetic field increases, it is clear that 
the interaction between thermally nucleated flux-line de- 
fects will be. less screened, driving the transition more 
continuoustll. The FLL melting transition thus becomes 
less first order in the anisotropic case as B increases, due 
to the presence of already thermally nucleated afe-plane 
vortex loops. Therefore the entropy jump associated with 
the transition is reduced as B increases. 

Recently, an interesting decrease in the entropy jump 
AS-jwith increasing B was in fact observed by Zeldov et 
alx3. Note that a consideration of fluctuations of the flux 
lines only, will give AS ~ B 1 / 2 , an opposite trendE-3. The 
above discussion shows that inclusion of vortex-loops into 
the consideration of the FLL melting in a natural way 
explains the data of Zeldov et al. The vortex- loop un- 
binding is crucial for the statistical mechanics of a type-II 
superconductor in magnetic fields B 6 (0 — 2T). 

We have discussed the existence of, and interplay be- 
tween, two distinct phase-transitions in type-II supercon- 
ductors. One is the melting of the FLL, which looses its 
meaning as the induction goes to zero. The other is a 



novel vortex-loop unbinding transition in finite magnetic 
field, with a zero-field counterpart. We emphasize two 
novel aspects of this work. The vortex-loop unbinding 
transition has been observed i) in a finite magnetic field 
and ii) both in the flux-line lattice and flux-line liquid 
cases, depending on our choice of T. 
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